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ABSTRACT 


Theoretical model and data analysis of DE-2 observations for determining 
the correlation between the neutral wave activity and plasma irregularities have 
been presented. The relationships between the observed structure of the 
sources, precipitation and Joule heating, and the fluctuations in neutral and 
plasma parameters are obtained by analyzing two measurements of neutral 
atmospheric wave activity and plasma irregularities by DE-2 during perigee 
passes at an altitude on the order of 300-350 km over the polar cap. Though it 
is still not clear whether one can conclude from these relationships that these 
perturbations arise from the Joule heating and precipitation, a theoretical model 
based on thermal nonlinearity (Joule heating) to give mode-mode coupling is 
developed to explore the role of neutral disturbance (winds and gravity waves) 
on the generation of plasma irregularities. 
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I. Introduction 

NASA awarded a research grant (Grant No. NASA - NAG5-479) bearing 
the title "A Morphological Study of Waves in the Thermosphere Using DE-2 
Observations”, to the Polytechnic Institute of New York with Professor Stanley 
H. Gross as the Principal Investigator for one year beginning September 1, 
1984. Under the support of this research grant, a number of investigations 
have been pursued. The investigation covered by the grant deals with data 
analysis of DE-2 observations for determining the correlation between the neu- 
tral wave activity and plasma irregularities. The aim of this study is to lead us 
to identify the source of perturbations, its structure and characteristics. More- 
over, it is also our intention to develop theoretical models for understanding 
the coupling between the neutral waves and the plasma irregularities and for 
interpreting wave properties obtained from measurements, so as to characterize 
and interpret the observed disturbances. 

At present, two measurements of neutral atmospheric wave activity and 
plasma irregularities by DE-2 during perigee passes at an altitude on the order 
of 300-350 Km over the polar cap are analyzed. They are for September 5, 
1981 at 21:13 UT (day 81248) and October 15, 1981 at 0:51 UT (day 81288). 
Magnetic activity on the first day was modest. The kP was 2 in the seventh 
three hour period (18:00 to 21:00 UT) and 3 + for the next three hour period 
(21:00 to 24:00 UT). The second case was at the start of a day that followed a 
day of considerable activity, October 14, 1981, in which kP reached 8. The kP 
index was 6_ from 21:00 to 24:00 UT on that day, but decreased to 3 from 0:00 
to 3:00 UT on October 15, 1981, which interval covered the time of the second 
case. The effect of the previous day's activity was therefore quite evident in 
the measured data. The calculated Joule heating flux was found to be about 
twice as great on the October 15 pass than it was on the September pass. 
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Of Interest here is the relationship between the observed structure of the 
sources, precipitation and Joule heating, and the fluctuations in the neutral and 
plasma parameters. Comparisons are made between perturbations of tempera- 
ture, densities, electric field turbulence, neutral winds, neutral densities, 
between wavelike behavior in neutral densities and irregularities in ion density, 
all as related to these source structures and as related to each other. Though 
some relationships are found, it is not clear, as yet, whether these perturbations 
arise from the Joule heating and precipitation. They may also be associated 
with the structure of the convection cells at the time, and measurements are 
needed to clarify this point. A preliminary theoretical model based on thermal 
nonlinearity ( Joule heating) is also developed to explore the role of neutral dis- 
trubance (winds and gravity waves) on the generation of plasma irregularities. 

In section II of this report, the background of the research area is dis- 
cussed. Section III presents the results of Data Analysis on the two measure- 
ments of DE-2 perigee passes. The theoretical study on the generation of 
plasma irregularities by collisional coupling of the neutral’s velocity perturba- 
tions to the ionospheric plasma is given in Section IV. 
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II. Background of the Research Area 
1.0 Background 

Wave-like perturbations of constituents in the thermosphere are detected 
by measurements of the ionosphere by ground facilities and by instrumentation 
on board satellites measuring both neutral and ionospheric parameters. Pertur- 
bations of ionization are known as traveling ionospheric disturbances (TID’s) 
(Hines, 1974; Georges, 1968; Yeh and Liu, 1974). A number of reports of 
large scale TID's (from hundreds to thousands of km) have suggested that they 
originate from high latitude regions during auroral substorms (Thome, 1968; 
Davis and de Rosa, 1969; Testud et al., 1975). These waves are believed to be 
caused by heat and momentum input into the atmosphere in the auroral zones 
from particle precipitation and Joule heating (Vickrey et al., 1982; Wickwar, et 
al., (1975). Their origins have been reviewed by Hunsucker (1982). Studies 
have also been made of these waves in association with magnetic storms by 
Richmond and Matsushita, (1975), Mayr and Volland, (1973) and Prolss, 
(1982). 

Medium scale TID’s (from tens to hundreds of km) are also observed, but 
their sources have not been well identified, though it is believed that meteoro- 
logical phenomena may produce many of them. Larsen et al., (1982) have 
associated thunderstorms with these waves. Mastrantonio et al., (1976) have 
also associated these waves with Jet streams in the atmosphere, and Bertin et 
al., (1978) have detected waves related to the Jet stream. Auroral sources are 
also believed to be the cause of many of these disturbances (Georges, 1968). 

The characteristics of all these are found to be generally in accordance with 
gravity wave theory (Hines, 1974) which has received some direct support from 
satellite measurements on board Explorer 32 (Dyson et al., 19700 and Atmos- 
pheric Explorer AE-C (Reber et al., 1975; Gross, 1980, 1984; Gross and 
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Huang, 1084; Gross et ah, 1080, 1081a and b, 1082, and 1084; Huang et ah, 
1080; Hoegy et ah, 1070). It is believed that the neutrals are perturbed by 
gravity waves as a result of localized sources of heating and momentum, and 
the neutrals, In turn, drive the ionization through collisions and temperature 
perturbations to produce these TID’s. Gravity wave theory is concerned with 
two types of waves called gravity waves and acoustic-gravity waves. Acoustic- 
gravity waves make up the upper frequency branch and gravity waves the lower 
frequency branch of the hydrodynamical system in which both gravity (buoy- 
ancy) and compressibility act as restoring forces. The upper branch becomes 
identical to ordinary sound waves in the audio frequency range, but at its 
lowest frequencies, such waves are called ’’infrasonic waves.” This branch has a 
lower cutoff frequency called the ’’acoustic cutoff frequency,” and there is a 
cutoff region below this frequency which separates the two branches. The 
lower end of the cutoff band, which is the upper limit of the gravity wave 
branch, is called the ’’buoyancy frequency.” The cutoff band is determined by 
the temperature and its vertical gradient and by the mean molecular mass of 
the medium, so that it varies somewhat with altitude. The cutoff band is rela- 
tively narrow, and periods within it range from about several minutes at low 
altitudes to from 10 to 15 minutes at higher altitudes, depending on tempera- 
ture. The periods of gravity waves then range from the period of the buoyancy 
frequency to many hours, whereas those of acoustic-gravity waves range from 
the acoustic cutoff period to periods of seconds and less. In the gravity wave 
branch, long period waves usually have large horizontal wavelengths and 
shorter periods smaller wavelengths (Hunsucker, 1982; Morgan and Tedd, 
1983). 

Though fluctuations in ionization in the thermosphere have been measured 
both on the ground and by satellite borne instrumentation, the detection of 
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wave-llke perturbations of neutral constituents in the thermosphere depends 
mostly on access by rockets and satellites. Radar is used only up to an altitude 
of about 100 km for neutrals. Satellite borne experiments have detected wave- 
like perturbations, both in total neutral density and temperature; for example, 
ESRO 4 (Prolss and Von Zahn, 1974a and b; Trinks and Mayr, 1976) OGO 6 
(Reber and Hedin, 1974; Taeusch et al., 1971), AEROS-A (Chandra and 
Spencer, 1975 and 1976; Trinks et al., 1976), AE-C (Reber et al., 1975; Potter 
et al., 1976; and the various papers authored or co-authored by Gross in the 
references) and DE-2 (Hoegy et al., 1981). The density variation of individual 
species, namely, nitrogen, oxygen, helium and argon, have also been measured 
on board ESRO 4 , AEROS-A, AE-C, AE-E, and DE-2. Vertical velocity vari- 
ations of these species have also been measured by AE-C, AE-E, and velocity 
variations by DE-2. 

2.0 Source Problem and Energy Distribution • 

2.1 High Altitude Sources 

Worldwide temperature increases in the thermosphere during magnetic dis- 
turbances have been observed by many researchers (Jacchia et al., 1967; Roe- 
mer, 1971; Reber and I-Iedin, 1974; Chandra and Spencer, 1976; and others). 
These disturbances are believed to be the source for many gravity waves 
observed in the thermosphere at mid and low-latitudes. The basis for this 
belief stems from observations of TID’s by ionosounders (for example, Klos- 
termeyer, 1969), by incoherent backscatter systems (Testud, 1970; Thome, 
1968; Testud et al., 1975) and by Doppler sounders (Georges, 1968) that detect 
these distrubances moving toward the equator with speeds ranging from 200- 
700 m/sec. The relationships between thermospheric neutral particle effects 
and magnetic activity is inferred as well from the correlations of measured data 
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with storm events (Trinks et al., 1075 and 1070; Chandra and Spencer, 1975 
and 1970). Waves have also been observed In association with auroral sub- 
storms in the evening section (Davis, 1971), though it has been suggested that 
some observed in the daytime may rise from the precipitation of soft particles 
about the cusp region (Trinks et al., 1975). 

Two mechanisms in the auroral zone are considered to be the sources for 
the waves. These are Joule heating due to electrojet current dissipation in the 
ionosphere and momentum disturbances arising from the Lorentz force (Chi- 
na on as and Hines, 1970; Cole, 1071; Testud, 1970; Richmond and Matsushita, 
1975; Roble et al., 1978; Richmond, 1079; Prolss, 1982) moving at supersonic 
speeds (Chimonas and Peltier, 1970) and heating due to particle precipitation 
(Hays et al., 1973). Nevertheless, it is generally accepted that the Important 
mechanisms at work are mainly Joule heating and the Lorentz force in the E 
region (see Hunsucker, 1982). Gross et al., (1984), however, suggest that 
heating by particle precipitation in the F-region may also be at work producing 
waves observed by AE-C and AE-E. The structure of diffuse aurora caused by 
lower energy particles (Lui et al., 1082) reinforces this idea. Furthermore, 
Hunsucker, (1983) finds that there is no apparent correlation between TID’s 
and magnetograms, except for large magnetic storms, suggesting that E region 
sources may not be at work when magnetic activity is below some threshold; 
New evidence of Morgan, (1983) may open such matters for further considera- 
tion. Anderson et al., (1982) propose that wind shears in the equatorial zone 
of the F region could also generate gravity waves in the thermosphere. 

The energy input to the thermosphere from high latitudinal regions is 
believed to be significant relative to solar EUV input to the thermosphere 
(Ching and Chiu, 1973). The high latitude input is therefore of great impor- 
tance to the physics of this layer, and the manner in which the energy is glo- 
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bally distributed to lower latitudes must be understood. Two major transport 
processes are favored; one Is based on the thermospheric winds, the other on 
waves. It Is very likely that both play a role (Klostermeyer, 1973; Chandra and 
Spencer, 1976; Mayr and Volland, .1972 and 1973). Waves distribute their 
energy through dissipation by viscosity, thermal conductivity, ion drag, wave- 
wave coupling and wind Interactions. The extent of transport of energy and 
momentum by waves is of interest here. Richmond and Roble, (1979) esti- 
mate 2-4 x 10 15 J can generate gravity waves capable of producing observed 
ionospheric disturbances at Millstone Hill and Arecibo incoherent scatter radar 
facilities. 

2.2 Low Altitude Sources 

Evidence that gravity waves in the thermosphere also originate in the 
lower atmosphere has been reported by Baker and Davies, (1969), Davies and 
Jones, (1971), Chimonas and Peltier, (1974) and Prasad et a!., (1975). These 
investigators have found connection between the thermospheric disturbances 
and severe thunderstorm activity at mid- altitudes. Others have also found 
correlations with ground level microbarograms (Bowman and Shrestha, 1966; 
Khan, 1970; Herron and Montes, 1970; Shrestha, 1971a and b). These latter 
meteorologically-connected events appear to have shorter periods than the large 
scale events, from about 3 minutes to the order of 1 minute. Waves of such 
periods are so-callecl acoustic-gravity waves of the upper frequency branch of 
the mode. Coupling to the stratosphere has also been noted (Mason, 1968 and 
1976; Fraser and Thorpe, 1976a and b). There have been observations of grav- 
ity waves in the thermosphere directly related to Jet stream activity (Gertin et 
al., 1978), confirming a suggestion by Mastrantonio et al., (1976). Though evi- 
dence of correlation of Jet stream activity and pressure fluctuations at the 
ground is plentiful, correlation with ionization variation would be an important 


- 9 - 


finding, because of implication as to the transport and distribution of energy 
and momentum. 

It may be expected that meteorological sources are relatively common, 
though the effects of an individual event may be localized (Davies and Jones, 
1971), in contrast with what is believed to be the global-like nature of distur- 
bances generated in the auroral region. However, the common presence of 
waves of scale sizes ranging from about 80 km to several hundred kilometers 
observed in satellite data (Gross et al., 1980 and 1981b, 1982; Gross and 
Huang, 1984) suggest that considerable meteorologically generated disturbances 
reach the thermosphere. In this regard, it is important to know what part of 
the disturbances are associated with auroral activity. Meterological events (not 
all of which appear to produce waves in the thermosphere) may cause the tran- 
sport of energy into the thermosphere in amounts not inappreciable relative to 
the solar EUV (Hines, 1965; Testud, 1970; Lindzen and Blake, 1970) If this 
much energy can reach the thermosphere, after possible partial reflection at 
some level and/or possible evanescence in some altitude region, such 
phenomena may also be indicative of the deposition of much more energy 
below the thermosphere. The chemistry and dynamics of the lower thermo- 
sphere, mesophere, and stratosphere can be significantly modified in such 
events. The interplay of tropospheric phenomena and these upper regions is of 
importance. One aspect may bear on the generation of turbulence and mixing 
up to the turbopause.. Wave effects on the turbopause may also be important 
because of the effect of the altitude of the turbopause on the specie distribution 
above It (Chandra and Herman, 1969; Chandra and Spencer, 1976). 
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III Measured Data and Analyses 

The orbital paths over the polar region are shown in Figure 1 for the two 
passes. These are plots of invariant latitude against magnetic local time. The 
path for the second clay for orbit 1071 on day 81288 passes closer to the mag- 
netic pole than that of the first day 81248, orbit 492. 

Figure 2 is a plot vs time in hours and minutes UT, as well as other 

\ 

ephemeris information, of the variation of various parameters measured by the 
LAPI, VEFI, LANG and WATS instruments for orbit number 492 on day 
81248. The top panel contains particle precipitation information in the form of 
the energy flux for electrons with energies as given by the logarithmic scale on 
the left. The flux is not readily obtained from the graph, since the original data 
utilized a color scale: however, the darker the appearance, the higher the flux. 
The highest possible flux of the scale is l erg/cm 2 - sec- sr- ev. The second 
panel is the turbulence in the electric field plotted on a logarithmic scale for the 
frequency band 8-16 HZ, shown as dots, and the 4-16 kHz band, shown as 0*s. 
The rise in the turbulence is quite evident at about 21:13 UT, at which time 
spikes in the precipitation flux of 3-4 kev electrons are quite evident in the 
LAPI data. These spikes are probably upward currents, and the region is 
believed to be the cusp. The third panel is a plot of electron temperature. The 
discontinuity and perturbation at 21:13 UT confirms the cusp location. The ion 
density is shown in the fourth panel, and a drop in density as well as irregulari- 
ties are observed to occur at about the same time. The fifth and bottom panel 
contains both the horizontal and vertical neutral winds, coded H and V, respec- 
tively. Note the rise in the horizontal wind at about 21:13 UT. Also note the 
wavelike variations in both wind components. The horizontal wind here is 
across the orbital pass and Is mostly east- west;. The step in the horizontal wind 
is not a true discontinuity. It is the result of the change in direction, as 
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designated, on passing over the geographic pole (note the geographic latitude 
80.98°). 

Figure 3 is the same as Figure 2, but for orbit 1071 on day 81288. The 
extended region of precipitaiton is quite evident and contrasts with that in Fig- 
ure 2. Note the spikey nature of the precipitation at about about 00;50 UT 
which is also associated with the rise in the electric field turbulence, the change 
in electron temperature and ion density there, as well as the rise in the horizon- 
tal wind. It is believed that this Is the cusp region. There is also considerable 
turbulence before this time which undoubtedly relates to the borad extent of 
the precipitation. The winds, electron temperature and ion density also exhibit 
considerable variation throughout, and the winds appear quite wavelike in their 
variation. 

Figure 4 contains plots of the variations in the magnetic field components 
and one electric field component, the x component which is along the space- 
craft velocity vector, vs time in hours and minutes UT for the two days, 81248 
on the left and 81288 on the right. Also on each of the graphs is a plot of the 
integrated downward Joule heating flux, designated SY in the figure. Note that 
the center of the flux is at 21:13 UT on day 81248 which is the time when the 
orbital path cuts through the cusp, and at 0:50 UT on day 81288, also the time 
of travel through that region. The peak Joule heating flux is about 50 mW/m 2 
on day 81248, whereas it is about double, 100 mW/m 2 on 81288. The duration 
of the Joule heating pulse in orbit at the 10% level, 5mW/m 2 , is about 53 
seconds, about 400 km in distance along the orbital track for day 81248. The 
Joule heating pulse for day 81288 is actually like a double pulse. Its overall 
duration at the 10% level, 12 mW/m 2 is about 120 seconds, or a distance of 
about 900 km along the orbital track. 

Figure 5 is a plot of the densities of oxygen, nitrogen, ions and electron 
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temperature vs time in seconds UT and other parameters. Samples are roughly 
1 second apart. [The upper graph is fo day 81248, whereas the lower is for day 
81288. These plots show the electron temperature and ion density variations in 
greater detail than those in Figures 2 and 3. For 81288 on the bottom, note 
the dip in the oxygen density, the rise in the nitrogen density, the increase in 
electron temperatures and decrease in ion density, all at about the same time, 
the time interval during Joule heating, about 0:50 UT. Note also the more gra- 
dual fluctuation of the neutral gas densities as contrasted with the irregularity of 
the ion density and electron temperature, the neutral density changes are not 
very evident for the Joule heating region 21:13 UT for 81248, though there is 
some small decrease in oxygen and increase in nitrogen densities. The more 
gradual fluctuation of the neutral densities in contrast with irregularities in the 
ion density and electron temperature is also quite apparent. 

Figure 6 is a plot vs time in seconds UT of the fluctuations in the densities 
of the ions, oxygen and nitrogen for day 81248. The fluctuations are obtained 
from the densities shown in Figure 5 by passing the logarithm of the densities 
through a digital bandpass filter with a lower cutoff frequency of 0.005 Hz, 
which corresponds to about 1500 km scale size along the orbital track, and an 
upper cutoff frequency of 0.833 Hz, which corresponds to about 90 km. Thus, 
the filter passes components of scale size between these limits. The upper 
cutoff is used to filter out noise in the data. Some of the features overlap. 
Thus, the negative excursion of the oxygen and electron density fluctuations 
and the rise in the nitrogen density fluctuation at a point past the 76350 second 
tick marker is associated with the peak in the Joule heating. Other both smaller 
and larger scale features match, though there are still some differences. 
Regions seem to correlate or anti-correlate in part. The fluctuations are the 
relative fluctuations with respect to the background. Relative ion density 
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fluctuations as large as -13% are obtained, though the mean level is closer to 
20% The oxygen-relative fluctuations are on the order of 3%, whereas that of 
nitrogen is closer to 3.5% 

Figure 7 shows the FFT power spectra of the fluctuations in Figure 6 plot- 
ted against wave period rather than frequency. The only spectral feature that 
appears common in the three spectra for this day, 81248, is one with a period 
of about 100 seconds which corresponds to a spatial structure of about 750 km. 
This size is roughly twice the scale size of the 10% level of the Joule heating 
pulse in Figure 4 for day 81248. This relationship could still be accidental. 

Figure 8 is the same as Figure 6, but for day 81288. The same filter has 
been used. The feature associated with Joule heating may be seen Just beyond 
the 2977 seconds tick mark. Other features appear to correlate to some extent, 
though there is still considerable differences. 

Figure 9 shows the power spectra plotted vs period for the fluctuations of 
Figure 8, which is for day 81288. Spectral features appear to correlate some- 
what at about 110 to about 128 seconds, 825 to 900 km on the orbital track. 
This scale size corresponds to the width at the 10% level of the Joule heating 
pulse. /In the following, we present an/ The connection may again be acciden- 
tal. The large peaks in the oxygen and nitrogen spectra at about a period of 
500 seconds is spurious and due to the finite size of the data samples. 


- 14 - 


IV Theoretical Analysis and Modeling Studies 

Modeling studies are necessary to guide the interpretation of measured 
data. Since the coupling of neutrals and ionization is of interest, it is necessary 
to take into account more than one fluid. Such studies have been made for the 
sourceless linear case (for example, Dudis and Reber, 1970; Gross and Eun, 
1976 and 1978; Del Genio et al., 1979). Since sources and their characteristics 
are of interest as well, source driven models must be treated. Such studies, 
however, were only made for atmospheres of Just one species (for example, 
Chlmonas and Hines, 1970; Richmond and Matsushita, 1975). At the 
Polytechnic we have been studying the problem of multiple constituents or 
fluids with sources, particularly for two species (Eun and Gross, 1976a and b; 
Gross and Eun, 1976, 1978). No extensive analytical treatment has as yet been 
made for more than two species, though with the assumption of truly minor 
species, one may treat the atmosphere as a two-fluid problem with one fluid for 
the major species and the other any of the minor species. Mayr and Volland, 
(1976) and Mayr et al., (1982 and 1983) have modeled the problem on a com- 
puter for propagation in a spherical geometry for a number of neutral consti- 
tuents with losses. No full treatment of ionospheric plasma in a magnetic field 
coupled to neutral species with sources has as yet been made. 

In the following, we present an analysis to show that plasma irregularities 
can be generated by collisional coupling of the velocity perturbations of the 
neutral waves in the neutrals to the ionosphere plasma. The thermal nonlinear- 
ity gives the mechanism for mode-mode coupling and provides a channel for 
transfer of the space-time oscillating perturbations into disturbances having only 
spatial oscillations (irregularities). In order to simplify the analysis, geometry 
corresponding to equatorial latitude is chosen. The application of the analysis is 
then focused on explaining the spread F (ESF) phenomenon. 
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1.0 Coupling of Nelural and Ionospheric Disb alances 

The rising and the subsequent descending of the ionosphere after sunset 
have been noted as the striking symptom of the onset of ESF (see, e.g., Her- 
man, 1066; Rastogi and Woodman, 1078). These upward and downward 
motions of the ionosphere are controlled by the E region dynamo. During the 
reversal of neutral winds, a wind shear may form and become the source of 
neutral waves (Hines, 1967, 1071). The formation of a wind shear in the F 
region was invoked by Anderson et al. (1082) as the hypothesized source of 
in-situ neutral waves for the local seeding of ESF. No wind shear in the F 
region was seen in the ALTAIR data, however, at time of wave structure for- 
mations. The wind shear required for the generation of in-situ neutral waves in 
the F region was estimated to be, at least, 18 m/sec/Km based on the criterion 
that R|< 0.25, where Rj is the Richardson number. While this required wind 
shear seems to be large in the F region, much greater wind shear can exist in 
the E region. 

Intense neutral waves produced, for instance, by the solar terminator 
(Beer, 1073b) or a wind shear in the E region are suggested in the present 
paper as the potential local coherent sources of exciting large-scale 
(> tens of kilometers) ionospheric density irregularities and forced ion acoustic 
modes* Ionospheric density irregularities are driven by the thermal pressure 
force that stems from the collisional dissipation of both the neutral wave and 
the excited ion acoustic modes via plasma-neutral collisions. This process, as 
shown below, relies on large pi asm a- neutral collisions. It thus represents an 

efficient coupling between the neutral and the ionospheric disturbances in the E 

< 

region rather than in the F region. However, the electric field perturbations 
associated with the large-scale field-aligned E region irregularities can map along 
the geomagnetic field lines up to the F region (Farley, 1060) and, consequently, 
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provide the initial plasma density perturbations for the subsequent excitation of 
Rayleigh- Taylor instability. 

1.1 Neutral wave as a pump 

Neutral waves, whose intensity is evaluated in terms of the neutral velocity 
perturbation (VJ in the following analysis, can transfer their momentum and 
energy to the charged particles through collisions. A plasma instability can be 
excited at the expense of the neutral wave. Whether the induced ionospheric 
disturbances can significantly affect the F region dynamics depends upon the 
characteristics of the instability such as the thresholds and the growth rates. 
The proposed neutral-ionosphere coupling process is analyzed by fluid equa- 
tions, because the scale lengths of plasma modes are greater than the ion 
gyroradii by several orders of magnitude. The nighttime E region is approxi- 
mated as a homogeneous plasma imposed by a uniform magnetic field. 

The two plasma modes under consideration are field-aligned low frequency 
mode and a forced ion acoustic mode that are parametrically excited by the 
neutral wave. This three-wave interaction process can be conveniently 
described by the following wave frequency (w) and wave vector (k) marching 
relations: 

R’e(k’o) = F^e( w s) ^e(^a) 

k = k + k 
^0 8 - — ’1 

where R e (o>) means the real part of u> the subscripts o, s, and a denote the 
neutral wave, the low frequency mode, and the ion acoustic mode, respectively. 
It will be shown that R e (cu s ) <<R e (a; a ) and the low frequency mode is nearly a 
zero frequency (or purely growing) mode. The forced ion acoustic mode is so 
termed because its frequency is determined by that of the neutral wave. There- 
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fore, this frequency may far deviate from the characteristic frequency 
(~a few kHz) of the ion acoustic waves In the ionosphere. 

The perturbations in plasma density (n), velocity (y), and electric field (E) 
caused by the low frequency mode (or ion acoustic mode) are represented by 
N s ( 6 N a ) , V ( 5VJ , and E ( <5 EJ /respectively. The plasma-neutral collision fre- 
quency has the expression of u pn — ^ pno (l + 5N n /N no ) with the neutral density 
perturbation (<5N n ) taken into account, where p = e (electrons) or i (ions), N no 
is the unperturbed neutral density, and u vrio is the collision frequency in the 
absence of <5N n . 


1.2 Large-scale Field- aligned Ionospheric Irregularities 

The coupled mode equation for the low frequency mode (i.e., the large- 
scale field-aligned ionospheric irregularities) can be derived from the following 
linearized electron and ion equations. 

They are: 

(I) Momentum equations 


M.N.^-%, = - V(T e0 N s + N 0 5T t ) + N 0 e(E, + X2 B 0 ) 


<5N, 


^e^o^enoYse ^•e^o eno 

^no 


eno 


( 1 ) 


MiN^V = - v(T te N, + N 0 «T,) + N 0 e(E + /v X2B 0 ) - 


«N. 


M N v V - M-N — 6V v + M- <5N*V v 

AVi i i>l o^moXsi AVA i AN o ~ai ^ino * iVA i vlN aI n i 

iN no 


ino 


( 2 ) 


where M e (Mj), e, T eo (T io ) , N 0 , and B (=2B C ) are the electron (ion) mass, the 
electric charge, the unperturbed electron (ion) temperature, the unperturbed 
plasma density, and the earth’s magnetic field taken to be the z axis of a 
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rectangular coordinate system; the last two terms on the RI-IS of (1) and (2) 
come from the coupling of the gravity wave (6N n and vj and the ion acoustic 
mode (^yj; the electron and the ion temperature perturbations (6T e and 6Tj) 
result from the thermal effect described by the heat source term in the energy 
equation; 

(II) Energy equations 

f N o {T « + N ° T *o V 4 V. 

M e 

= V(R||V||+Rx c Vx)6T e -3N 0 P en0 — (<5T e -6T n )-n 0 i/ en0 M e («V£- V) 

(3) 


_N — <5 T 4* N T X7 • V 
2 0 c>t i ^ ° i° V 

Mi 

= V *(R|V v ||+RxV x) «T r 3N 0 P in0 - (<$T r 5TJ- N 0 i/ ln0 M, ( 6 V* • V ) 

(4) 

where M n and 6T n are the neutral mass and the neutral temperature perturba- 
tion, respectively; the first term on the RHS of (3) and (4) represents the 
parallel and the cross-field heat conduction rates, the second term is the cooling 
rate due to plasma- neutral collisions, and the last one is the heat source term 
contributed form the collisional dissipation of both the gravity wave and the ion 
acoustic mode; 


(III) Continuity equations 

“N s + V ’ N 0 V = 0 = -?- N s +V • N 0 V (5) 

dt ■ ~ se dt 

where the quasi-neutrality condition, viz., N se c^N si = N s , has been used; 
further, it is reasonable to assume that 
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(IV) 


6TjC^6T n 


(0) 


if the ion mass is not very different from the neutral mass. Therefore, the 
cooling term in the ion energy equation disappears with this assumption. 

We also assume that the perturbations associated with the low frequency 
mode have a space-time dependence of the exp[i(k s x-w s t)] form, where 
w s = w r + i7 and is the real frequency (growth rate) of the low frequency 

mode. This mode has a lie Id- aligned nature, viz., its wave vector is orthogonal 
to the earth’s magnetic Held (B = 2B 0 ) and is chosen to be the x axis of the 
coordinate system. Replacing v( = 3cd/dx) and d/St in (l)-(6) by ik( = ixk s ) 
and - iw s and eliminating <5T e ,<5Tj , Y se ,V si ,andE from these equations, we 
obtain the following coupled mode equation for the field-aligned ionospheric 
irregularities 


( M e C 2 g 

(7— lw r )['Y-iWr + U ino + — — U eao ] + Cl j( 
l Mj 1 + V. 


2 K , T e 


+ 3 Mj ^ 7 


U + 7i 
M. -5V* 


M // T 


L n 7e 7i 
M a U onn 


1 2 F f \ . fe* ^ ^ino c * *] 

k S 2 ^eno(TT-) — + (2— ~ + 1)— *Yai * Y n 

L Mj 7 e M n 7 e 7j J 
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where V te (V ti ) and u e (u { ) are the electron (Ion) thermal velocity and the 
effective electron (ion) collision frequency, respectively. The three terms, 

2(M e /M n )^ ctl0 , ^ e k s 2 V‘/(n e 2 + u;), and v x k 2 i^/(n 4 2 + i/*) 

are the electron cooling rate due to electron-neutral collisions, the electron 
cross-field heat conduction loss rate, and the ion cross-field heat conduction 
loss rate, respectively. The RHS of (7) represents the nonlinearity of causing 
ionospheric irregularities, viz., the collisional dissipation of the gravity pump 
wave and the forced ion acoustic mode. 

1.3 Forced Ion Acoustic Mode 

The forced ion acoustic mode is the "high frequency” sideband excited 
parametrically by the gravity pump wave. Since its frequency is determined by 
the beat frequency of the gravity wave and the low frequency mode, it is not 
necessarily the characteristic frequency of the acoustic waves. The electron and 
ion equations that are required for the derivation of the coupled mode equation 
for the forced ion acoustic mode include the continuity equations and the 
momentum equations. 

For simplicity, the forced ion acoustic mode is assumed to propagate along 
the earth’s magnetic field; viz., = zk a . This assumption is quite reasonable in 
the E region where u it ^ 0 >>Cl i indicating that the ion dynamics can be treated 
approximately as in the. unmagnetized case. With the exp[i(k a z - cut)] type of 
perturbations, the linearized electron and ion equations can be written as: 

(I) Continuity equations 

- cu a <5N a + iN 0 k a 6V aez = 0 =- icu a <5N a + iN 0 k a 5V aiz (8) 

namely, 
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*V MI = < 5 V aiz = < 5 V az , 


(II) Momentum equations 


M e N 0 ( - iw + P eno ) 6V a 


ae 


= IZk a Teo 5 N a - ' N 0 e(Z 5 E a + — ^ x ZB 0 ) 
6 N n 

V 4- M NV V 

N ~se ~ lvl e 1N s ^eno-Zn 


^e^o^eno 


•no 


M.No^lw + O < 5 V ai 


( 9 ) 


iZ a T i0 < 5 N a + N 0 e(Z6E a + - < 5 V ai x ZB 0 ) - 


MjNoV, 


< 5 N 


n 


mo 


V* + M-N !/• V 

tvt -si T iV1 » iN s ^ino^n 
iN no 


(10) 


where u> = iu a -f i7 . The quasi-neutrality approximation, 6Nj~ 6N e = 6N a , 
has been used in (8) - (10). 

Eliminating V sez , V siz , and 6E a from the z components of (1), ( 2 ), ( 9 ) 
and (10) leads to 


(“ “ ( w a _lc a C s") + (J'ino + T7" ^eno) “ 
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2 ■ 


7 + p. 


ino M e P eno 

T — 


ino 


6 N„ 


N, 


no 


M e 

^^az (^ino ^ ^eno) 


n; 

N, 


V, 


nz 


Pi 


ino 
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l iv 


■^e^ino 


7 "1* ^ino ^i(7 “I" ^eno) 


eno) J 


< 5 N a 


N 0 

V n V 


(ID 


^Nn 

N„„ 


with the aid of (8), where C s is the ion acoustic velocity. Under the assumi>- 
tion that R c (a; a i >7 and 7«P ino <P eno and since 


^a- KC* 


then (11) reduces to 



6V 0 


N£ 

N„ 


V 


nz 


1 - 


6N* 


k V* 
"a v nz 


N, 


cu. 


- 1 


NsV 

N„ 


nz 


( 12 ) 


Solving (9) and (10) for , 6V MX .^V^y and 6V ajy , respectively, we get 

the following results: 


u pn 0 • N,* u pnr , N,* 

*V Mx ^(^) 2 (-^V nx )- ( ““) ( TT"V ny ) 


n 


N, 


n e N 0 


(13) 


v N' 

^eno N o / 1 s 


^eno . , N s 


,5V ~ (— — ) 2 (—— V ) + ( — — -1 ( — — V ) 

u v aey — \ ^ ^ N ny ' ^ ^ v nxJ 


n 


N r 


(14) 


N s * N s * u eno 

6V aix ~ a — V nx + 0 — i_(_i]12.)v 
N 0 N c 'fie 


ny 


(15) 


and 


N s * 

Sv • / ' v - y Oi 

u v aiy — « 


N s * 

V„ v - fi — v, 


N, 


nx 


(16) 


where a = i and 0 = fi j/j/ ino for £1 j« ^j n0 in the E region; a = (v ino /Cl j) 2 
and 0 = t/ ino /f 1 i for Cl j»i/ ino in the F region. Since (i^ en0 /Cl e )~10 -2 (1CT 6 ) 
and (Cl i /J'j no )~ 10~ 1 (l° +2 ) in the E(F) region, Equations (12)-(16) yield a 
simple coupled mode equation for the forced ion acoustic mode, viz., 


<5V. 

— ai 


N* 

iN s 


N n 


y„ » 


(17a) 


for Q i « v. no in the E region, or, 


^Yai-^YniXz) Ko/^i) 


N s * 


N, 


+ Z 


N* 


N„ 


V„z »5Yae 


(17b) 


for fl j » u [no in the F region. It has seen that if the neutral wind blows dom- 
inantly along the earth’s magnetic field (i.e., V ~ ZV nz ), 6V ( in the F region is 
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2.0 Excitation of the instability 

Under the assumptions that u eao > u ino » 7 and k 2 C s 2 » cu 2 , can be 
separated into the following two equations 


V 


U) T ~ 


nx 




r OalAs) (KlAa) 


R e ( w a) « R e (w a ) 


(18) 


and 
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eno 
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M i 7e 
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~n 


(19) 


According to (18) and assuming that R e (cu a )> 7 , cu r (=R e (u; s )) is less than 
R e (w a ) by at least four orders of magnitude since 


w al/ lc s~ Kl/ k a~ c s~10 3m /s e c. 


Therefore, this result shows that the low frequency field-aligned mode can be 
considered to be a zero-frequency mode and that the frequency (R e (^a)) °f the 
forced ion acoustic mode is determined by that (cu 0 ) of the neutral wave, i.e., 

^•e( w a) ^* 0 * 

Based on (I7a,b) and the assumption that v ino / 7 f i and v eno /^f e »l, It 
becomes clear in (19) that the neutral- electron interaction ( *^n) k 33 

a negligibly small effect on the excitation of the instability as compared with the 
neutral-ion interaction ( ?za/ ino 5y^ *V). Hence, the RHS of (19) can be approxi- 
mated as 


v 


1 kH2 ^^. + d zam-sv- .v 
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showing that the proposed mechanism depends on the effective Ion-neutral 
coupling. Substituting (17a) Into (19) leads to the dispersion relation of the 
Instability, which Is a cubic equation of *7 (the growth rate), viz.. 

An 3 + B 7 2 + 4- D = 0 ( 20 ) 

where 

-A. 1 + t , eno t , j no /fl e fl j, 

B = 2(M e /M n )t' eno + VjK*Vll{n ?+V?) + e+t' 2 ) 

+ (io/3)i/ eno k s 2 v u Vn e n it 

c = ls iao V eao [k 2 V 2 /(ft j 2 +t'j 2 ) j(2M e /M n + ks 2 v£/n e 2 ) 

+ (2/3)^ cno (ks e n j) jlO(M e /M n )i/ eno 

+ «. 2 + 4i/ lno k s 2 vJ/(n , 2 +^)] - <2/3)(k 5 s <w'mo/n .n i) K I 2 , 

and 

D = - (2/3)i/ ino i/«; 0 (k=/n ,n j)( 4 M e /M n + k 2 v 2 /a *) K I 2 + 

2^too^,lo(k 5 XVn ,n ,) [kfvj/fn , 2 + ^ 2 )](2M,/M 11 +k 2 v 2 /n 2 ) . 

Although ( 20 ) is the dispersion relation derived for the E region case, it can be 
also applied to the F region when V~2V nz , namely, when the neutral wind 
blows primarily along the earth’s magnetic field lines. 


2.1 The threshold condition 

The threshold condition of the instability that is obtained by setting 7=0 
in the dispersion relation is given by 



3^ f M e k/V 2 ] f M e k^V 2 1 

(Bf + Vf) M n fl 2 j + n 2 


( 21 ) 
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where the two terms, 2M C /M n and k s 2 V t 2 /n 2 , correspond to the electron cool- 
ing due to electron-neutral collisions and the electron cross-field diffusion, 
respectively. The dormer effect dominates over the latter when the scale 
lengths (X s = 27rk s _1 ) of field-aligned ionospheric irregularities exceed 
(27rV te /n e ) (M n /2M e ) 1/2 ~ 8 meters in the ionospheric E region. 


Since X s (=27rk s 1 ) is of the order of the gravity wavelength 
(^.tens of kilometers), (21) can be approximately written as 


hf ,0 ,3, (27rV$)- 

l V th I'— ( ~) 7—7 ^ 

2 (fi +i/ f )\ s 


27TV 2 


(f-) ( v 

2 ^ino^s 


) ! 


3, , 27rV tU 8 


(— ) ( 
K 2 JK O.X 


■r 


r's 


(E region) 

(F region, when 

V ~zV 1 
~n — L v nz/ • 


This expression shows that the threshold condition of the proposed instability is 
primarily determined by the ion cross-field diffusion damping, whose exact 
form has been retained. In the ionospheric E region, the ion effective collision 
frequency is dominantly contributed from the ion-neutral collisions (p in0 ) 
that are much greater than the ion gyro-frequency (O j), that is, v\~v ia0 »n j. 
This means that ions are essentially unmagnetized and their cross-field diffusion 
damping is thus identical to their parallel diffusion damping. If we use the fol- 
lowing E region parameters (Rishbeth and Garriott, 1969): 


M(NO + )/M e = 5.5 X 10 4 ,R i/27T = 25.4 Hz, ^ ino (^eno) 

= 5.8 X10 3 ( 9.2 X 10 4 ) Hz, V ti = 2.54 X 10 2 m/sec 

(or Tj=210° K) and X s = lOOKni, then |V th |~ 8.5 x 10 _4 m/sec. By contrast, 
fVtjjl is independent of u ino in the F region because Cl s » £> ino in the case that 
\^~zV nz . This indicates that ions are strongly magnetized and their parallel 
diffusion damping may greatly exceed their cross-field diffusion damping if the 
irregularities have a non-zero parallel wave number. With the substitution of 
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the typical F region parameters: 

M(0 + )/M e = 2.9xl0 4 ,B j/27T = 4.76Hz,V u = 7.6xl0 2 m/sec(Tj = 1000°K), 

and X s =100Km, the threshold neutral velocity perturbation is 1.5 x 
10“ 1 m/sec that is greater than IV^ | in the E region by about two orders of 
magnitude. But, if the neutral wind blows dominantly across the earth’s mag- 
netic field (i.e., \^ t it can be expected from (17b) that the threshold velo- 

city in the F region is enhanced by another factor of B j/i/ jno ~400. In other 
words, 60m /sec) in the F region exceeds that in the E region by four 

orders of magnitude. The threshold in the E region seems to be rather small 
compared with the typical neutral velocity perturbation of gravity waves that is 
of the order of 40 m/sec. However, a small threshold does not necessarily 
ensure the growth rate to be large enough for the practical excitation of the ins- 
tability. As shown later, the proposed instability can be only practically trig- 
gered by gravity waves with periods of a few tens of minutes. 

2.2 The growth rates 

For the excitation of modes with X s »( 27 r V^/fi e ) (M n /2M e ) 1/,2 ~ 8 (13) 
meters in the E (F) region, the coefficients of the cubic equation (20) have the 
following simple form: 

Aosil + ^eno^ino/^ i (E region ) and 1.0 (F region). 

Bc=:2 (M e /M n ) P eno ~ 3.3 (E region) and 4.1xl0“ 4 (F region), 

C^2^ eno (k s 2 V 2 /^ ino ) (M e /M n ) [1- (l/3)(i/ ino /n i) 2 (V n /V ti ) 2 ] 

'-'-'1.42x10“ 5 (E region) and l.Ox 10 -6 (F region), 

D~ 4(k s V te i/ eno /n e ) 2 (k s 2 V 2 /p ino )(M e /M n )[l- (Vh/V^) 2 ] 

~ - 1.48x l0“ 5 (E region) and - 1.86x 10" 12 (F region), 

for V n = 40m /sec, X s = 100 Km, ^ ino =5.8 x 10 3 Hz(0.5Hz), and u tVl0 = 9.2 xlO 4 
Hz (12 Hz), in the E (F) region. Therefore, instead of solving the cubic 
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equation, 7 can be derived approximately from 
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where V th is given by (22). Since V n »V th , Equation (23) can be written as 
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for X s = 100 Km. It is apparent that the neutral-plasma interaction is much 
more effective in the E region than in the F region. The growth rate is 
Inversely proportional to the scale length of field-aligned ionospheric irregulari- 
ties that is of the order of the neutral wave wavelength. While a maximum 
scale length of 400 Km was detected by ALTAIR from the horizontal wave 
structure, the nighttime ionosphere spatially modulated with shorter scale 
lengths (~35 Km) prior to the onset of ESF was sensed by the amplitude scin- 
tillations of gesostationary satellite signals at 136 (MHz) (J.A. Klobuchar, 
private communication, (1984). The growth rate in the E region for 
X s = 35Km is 6.0 x 10" 3 sec -1 in this situation. Since the wavelengths of neu- 
tral waves are tens of kilometers or longer, the instability is expected to have 
10 - 2 sec -1 as the upper bound of its growth rate in the ionospheric E region. 
The growth rates have to exceed the molecular ion recombination rate 
(~ 10 -3 sec -1 ) for the operation of the instability. 

In the formulation of the theory, R e (w a )j >.7 has been assumed, where o> a 
is the angular frequency of the forced ion acoustic mode. The real angular fre- 
quency, R e (u> a ), is essentially equal to the neutral wave frequency (w 0 ) because 
of the excitation of ionospheric irregularities as nearly zero-frequency modes. 
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This condition requires that the neutral wave frequency (cu 0 ) be greater than 
7~2.0 x 10“ 3 1-Iz, namely, the neutral waves have periods (T 0 = 2n /u> 0 ) less 
than one hour. A clear portrait of the neutral pump waves can be delineated at 
this point: they are characterized with short periods (< one hour) and 
wavelengths of tens to hundreds of kilometers. 

0.30 Seeding the ESF 

The neutral wave-plasma interaction, as shown above, is not effective 
enough for the excitation of large-scale ionospheric irregularities in the F 
region. One may note, however, from (24) that in the F region, 7 oc i/ " * /2 and 
argue that smaller P ino at higher altitudes can lead to a larger growth rate. Two 
other things should also be noted: (1) p jno cannot be too small because 
p ino »7 has been assumed in the formulation of the theory, and (2) 7 ec v 
and P eno is smaller at higher altitudes. Actually, P eno /P ino < 10 2 at any altitude 
in the ionosphere. Hence, even if we talce i'enoA'ino = 10 2 , the growth rate is 
5.2x10“ 4 sec -1 for \ s = 35 Km in the case that the neutral wind blows along 
the earth’s magnetic field. This growth rate can exceed the atomic ion recom- 
bination rate above a certain height. Nevertheless, if the neutral wind blows 
dominantly across the earth’s magnetic field, the growth rate in the F ; region is 
reduced by at least two orders of magnitude. Moreover, the conditions for the 
instability become extremely stringent when the parallel heat conduction loss in 
the F region is not negligible. This effect can largely raise the threshold and 
decrease the growth rate of the instability. The excitation of the instability with 
a rather small growth rate requires almost a continuous source of neutral 
waves. 

While the forced ion acoustic modes are highly damped in the E region, 
the concurrently excited large-scale ionospheric irregularities have much less 
diffusion damping and can survive for hours after the disappearance of the 
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neutral-wave source(s). The electric field perturbations associated with the 
large-scale E region irregularities can extend along the geomagnetic field and 
many onto the F region (Farley 1060). These intense electric field perturba- 
tions then seed the F region for the subsequent excitation of the Rayleigh- 
Taylor instability that has been generally accepted to be the cause for the spread 
F echoes on the ionograms and the plumes on the VHF backscatter radar maps. 
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IV. Summary 

The data measured by the Dynamics Explorer 2 spacecraft (DE-2) provide 
a unique opportunity to determine many geographical parameters. In this 
report we present the results of analyses on two sets of data measured in two 
different days. On comparing the data, it is seen that there is a lasting effect of 
previous magnetic activity on the fluctuations for day 81288 in that the fluctua- 
tions are about double those for day 81248. This doubling is consistent with 
the Joule heating for the two days. The power spectra are also in agreement 
with this ratio. The ratio of powers being the square of the fluctuation ampli- 
tude is found to be about 5.7 for oxygen and 4.0 for nitrogen. Our analyses 
appear to favor Joule heating as the driver of observed wavelike fluctuations in 
neutrals rather than precipitation. It may be that the neutrals are excited by the 
entire auroral oval in the Joule heating pulse and filtering yields wavelike, 
rather than irregular, variation. Ion waves are then excited by the neutral 
waves through the ion-neutral collisional coupling. 

A theoretical model exploring the role of neutral waves on the generation 
of large scale plasma density irregularities is also presented. This model is prel- 
iminary in the sense that equational latitude is considered so that the effect of 
anistropy introduced by the geomagnetic field on plasma motion is minimized. 
This simplifies the analysis a great deal. .Nevertheless, the results of the 
analysis can be applied to explain the the equatorial spread F (ESF) 
phenomenon. The neutral waves of interest are characterized by wavelengths 
of tens to a few hundreds of kilometers and by periods of a few tens of 
minutes. The efficient excitation of ionospheric irregularities by our proposed 
mechanism relies on large plasma-neutral collisions. Such an optimum 

r . . ■ 

environment can be provided by the E region rather than the F region. Since 
neutral waves interact primarily with ions, the threshold of the instability is 
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determined by the ion heat conduction loss. It is shown that the plasma irregu- 
larities can be generated within a few minutes when large amplitude gravity 


waves are present. 
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Fig. 1 


Fig. 2 


Fig. 3 


Polar plots of the orbital passes for September 5, 1981, 
day 81248, orbit number 492 and for October 15, 1981, 
day 81288, orbit number 1071. The plot is invariant lati- 
tude vs magnetic local time. 

Plots of electron precipitation (toppanel), electric field tur- 
bulence (second panel) 8-16 Hz as dots and 4-16 kHz as 
■' the letter O, electron temperature (middle panel), ion 
density (fourth panel) and horizontal (H) and vertical (V) 
winds vs time in hours and minutes UT, invariant lati- 
tude, magnetic local time, altitude, local solar time, lati- 
tude and longitude. This figure is for orbit 492. on day. 
81248 for the twelve minutes from 21:10 to 21:22. 


Plots of electron precipitation (top panel), electric field 
turbulence (second panel) 8-16 HZ as dots and 4-16 kHz 
as the letter 0, electron temperature (middle panel), ion 
density (fourth panel) and horizontal (H) and vertical (V) 
winds vs time in hours and minutes UT, invariant lati- 
tude, magnetic local time, altitude local solar time, lati- 
tude and longitude. This figure is for orbit 1071 on day 
81288 for the twelve minutes from 21:10 to 21:22. 
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Fig. 4 


Fig. 5 


Fig. 6 


Graphs for day 81248 on the left and day 81288 on the 
right. Each graph contains plots of the measured variation 
in the x, y, z components of the magnetic held relative to 
a Magsat model ABX, ABY , ABZ , the x component of 
the electric field EX and the integrated downward Joule 
heating flux SY vs time in hours and minutes UT for each 
day. x is along the velocity vector of the spacecraft, y is 
upward and z is eastward or westward according to the 
spacecraft direction, north or south. The plot SY has 
been plotted ,as positive to mean downward, opposite that 
of the y axis. 


Graphs of the densities of oxygen, nitrogen and ions and 
the electron temperature vs. time in seconds UT as well as 
hours, minutes and seconds UT, shown above each graph. 
Also given are scales for altitude, latitude, longitude, solar 
zenith angle, local magnetic time and invariant latitude in 
each graph. The density scales in particules/cm 3 are loga- 
rithmic and on the left and the temperature scale in °K is 
linear and on the right. The upper graph is for 81248. 
The lower for 81288. 

Plots of fluctuations obtained by passing the logarithm of 
the ion, oxygen and nitrogen densities through a bandpass 
filter that passes structures with scale sizes from 90 to 
1500 km. The plots are vs time in seconds UT. The data 
are for orbit 492 on day 81248. 



Fig. 7 


Fig. 8 


Fig. 9 
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FFT spectra of the fluctuations in Figure 6 for day 81248. 
The spectra are plotted vs period in seconds. 


Plots of fluctuations obtained by passing the logarithm of 
the ion, oxygen and nitrogen densities through a bandpass 
Alter that passes structures with scale sizes from 90 to 
1500 km. The plots are vs time in seconds UT. The data 
are for orbit 492 on day 81288. 


Plots of fluctuations obtained by passing the logarithm of 
the ion, oxygen and nitrogen densities through a bandpass 
Alter that passes structures with scale sizes from 90 to 
1500 km. The plots are vs time in seconds UT. The data 


are for orbit 492 on day 81288. 



Fig. 1 



Polar plots of the orbital passes for September 5, 1981, 
day 81248, orbit number 492 and for October 15, 1981, 
day 81288, orbit number 1071. The plot is invariant lati- 
tude vs magnetic local time. 
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Plots of electron precipitation (top panel), electric field 
turbulence (second panel) 8-16 Hz as dots and 4-16 kHz 
as the letter 0, electron temperature (middle panel), ion 
density (fourth panel) and horizontal (H) and vertical (V) 
winds vs time in hours and minutes UT, invariant latitude, 
magnetic local time, altitude, local solar time, latitude and 
longitude. This figure is for orbit 492 on day 81248 for the 
twelve minutes from 21; 10 to 21; 22. 
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Fig. 4. Graphs for day 81248 on the left and day 81288 on the right. Each graph contains 
plots of the measured variation in the x, y, z components of the magnetic field • 
relative to a Magsat model ABX, ABY, ABZ, the x component of the electric 
field EX and the integrated downward Joule heating flux SY vs time in hours and 
minutes UT for each day. x is along the velocity vector of the spacecraft, y is 
upward and z is eastward or westward according to the spacecraft direction, 
north or south. The plot SY has been plotted as positive to mean downward, 
opposite that of the y axis . 
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Fig. 5. Graphs of the densities of oxygen, nitrogen and 
ions and the electron temperature vs time in 
seconds UT as well as hours, minutes and seconds 
UT, shown above each graph. Also given are scales 
for altitude, latitude, longitude, solar zenith angle, 
local magnetic time and invariant latitude in each 
graph. The density scales in particles/cm^ are 
logarithmic and on the left and the temperature 
scale in °K is linear and on the right. The upper 
graph is for 81248. The lower for 81288. 
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• Plots of fluctuations obtained by 
passing the logarithm of the ion, 
oxygen and nitrogen densities through 

* a bandpass filter that passes struc- 
tures with scale sizes from 90 to 
1500 km. The plots are vs time in 
seconds UT. The data are for orbit 
492 on day 81248. 


Fig. 7. FFT spectra of the fluctuations in 

Figure 6 for day 81248. The spectra 
are plotted vs period in seconds. 



















